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SUMMARY 

The total dose in rads-tissue from solar protons was tabulated for weekly time intervals, and the 
number of weeks which gave a dose above 25 rads behind 10 g/cm 2 of aluminum for the active 6 years of the 
19th cycle were called dangerous or large event weeks. The number of such event weeks was found to be 
only 3 weeks for the past 20 years. Even though the chance for smaller events is examined, it was found 
that for any reasonable high confidence level (95%), the smaller events could be ignored. Consequently, 
the total particle flux for the 19th cycle was divided by a factor of 3 and determined a single large 
event week. Using this spectrum, the tissue dose in rads is calculated at the center of an aluminum 
spherical shell. To correct for geometric effects and self shielding, this dose may be reduced by a 
factor of about 3. To predict the probability of an event occurring, the Poisson distribution was the 
most logical choice. The confidence one can use in employing the Poisson process and arriving at confi- 
dence levels for the experimental value of the mean is investigated. Several examples are given for dif- 
ferent mission lengths, and comparisons are made to other authors’ results. An extension of the Poisson 
process is made to incorporate the concept of small sample theory and arrive at the expected distribution 


function which answers the following question. If x Q 
bility of seeing x events in any observation time t? 

INTRODUCTION 

This study is concerned with the practical 
treatment of hazards from solar proton events out- 
side the magnetosphere of the earth. It is not con- 
cerned with prediction of flares as such nor the 
long-range solar cycle indicators such as sun 
spots. The study of all solar proton events as 
such is only casually related to the problem of 
large dose rates inside realistic spacecraft. For 
example, in reference 1, there are 76 events listed 
from 1942 through 1963. Since 1963, there are pro- 
bably an additional 24 more yielding a statistical 
sample of some 100 events. However, the corre- 
lation of this large sample with the events which 
are of real danger to space flight is very poor. 

For example, only 6 large events from 1950 to 1969 
would have given about 85% of the total 20-year 
proton dose behind a thin wall of only 2 cm of alu- 
minum. In addition, these 6 flares occurred in 3 
weekly periods . Thus . 3 of these 6 large events 
occurred from July 10-16, 1959, 2 occurred from 
November 12-18, 1960, and on February 23, 1956, 1 
large event was observed. The duration of an event 
is from 1 to 3 days. Even with this small sample, 
it seems that if conditions are suitable for a 
large event, the odds are very good that it will be 
followed within hours by another large event. Per- 
haps a time span of at least 1 week should be used 
to depict a total solar event, and it would be des- 
ignated as the solar proton flux or dose per week. 
Under this definition, there are only 3 sample 
weeks of large solar events from about 1950 through 
1969 (1,000 weeks). We are not dealing with an 
ordinary problem of statistical analysis, but with 
rare large events. 

The reason that the writer has undertaken the 
awesome task of predicting the improbable is not by 
virtue of his background in solar physics or sta- 
tistical analysis but because of his concern with 
protection of man and his radiation-sensitive 
equipment from space radiations for realistic 
spacecraft and missions. However, the reader should 
ask, why not leave this field of statistical astro- 
nomy to the experts. The answer is that the en- 
vironmental scientists do not have to design or 
evaluate realistic shields for sensitive film or 
radiation-conscious astronauts. 


events are observed in time t 0 , what is the proba- 


The consequence of having several solar pro- 
ton prediction models (which I do not wish to 
evaluate) has led to a wide disparity in results, 
especially when a reasonably high level of statis- 
tical certainty is desired. For example, at the 
99% probability (percentile) level, the predicted 
dose behind 20 g/cm 2 of aluminum for a l~year 
mission may vary by a factor of 10 or more between 
different writers Now this may not sound too 
bad considering the nature of the problem, but if 
a mission is planned with the requirement that the 
astronauts should not receive doses exceeding 100 
rads skin dose with a 99% probability, and if one 
solar proton dose prediction model requires a 
shield of only 15 g/cm 2 and another model 50 g/cm 2 , 
the radiation analyst has to make a vital decision, 
possibly affecting the life of an astronaut. Since 
no one wishes to be responsible for making a de- 
cision which could lead to dire consequences, the 
most pessimistic model is often chosen. This 
writer would not question this approach except 
that the desire to be on the safe side may readily 
get out of hand with the subsequent loss in mission 
capability. 

A PRIMITIVE DOSE MODEL IS PROPOSED 

The most direct measure of the hazard of a 
given solar proton event is the rads (.01 J/Kg) - 
tissue absorbed dose that would be measured behind 
various thicknesses of a typical spacecraft mater- 
ial such as aluminum. The simplest method is to 
find the point tissue dose at the center of a 
spherical shell of aluminum. However, due to self 
shielding by the astronaut (approximately a factor 
of 2) and spacecraft geometry as well as on-board 
equipment, the actual skin dose may be less by a 
factor of 3 or more than the point dose at the 
center of a spherical shell. This factor varies 
depending on the solar proton spectrum, the space- 
craft geometry, and the location of the astronaut 
in the craft. This writer suggests using a factor 
of about 3 reduction for the solar proton point 
dose at the center of a spherical shell in order to 
estimate the likely skin dose to an astronaut in a 
real spacecraft of a given average thickness. (It 
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should be noted that this factor of 3 may be too 
large for a very "hard" spectrum such as seen in 
the trapped radiation belts, when a factor of 2.5 
may be more in order.) To correct for depth dose 
(bone-marrow depth) , a thickness of 5 cm of tissue 
is often employed. This is approximately equal to 
6.5 g/cm 2 of aluminum in equivalent shielding 
effectiveness . 

In order to clarify the exact assumptions 
which are to be employed, the following infor- 
mation is pertinent. The only adequate data avail- 
able at this writing for solar proton predictions 
are from the 19th cycle (1954-64). However, there 
have been observations of sun spots (indicator of 
solar activity) for about 200 years. Based on 200+ 
years of observations, where the first cycle would 
date back to the middle of the 18th century 
(average cycle length is about 11 years) , the 19th 
cycle had the highest maximum sun spot count yet 
recorded. In fact the average maximum sun spot 
count is more than a factor of 2 lower. The 20th 
cycle, which we are now well into (past the peak 
activity) , has a sun spot count somewhat above this 
average. There has not been a large solar proton 
event comparable to the eight largest events of 
the 19th cycle (dose behind 10 g/cm 2 ). At this 
point, one may be 'led to believe that the 19th 
cycle is a fairly* rare type cycle. With the pre- 
sent low occurrence of large solar proton events 
and the rather extensive sun spot counting dating 
back over 200 years, one might conclude that the 
probability of getting a solar cycle as active as 
the 19th is on the order of 1/20 or 0.05. This, 
of course, cannot be objectively demonstrated, and 
will not be, unless considerably more knowledge is 
obtained about the physics of the sun. One valid 
objection to the above is that the sun spot number 
is a poor indicator of large solar proton events. 
Also, the sun spot indicator may have changed 
during the last 50 years due to better observa- 
tions, so that possibly the first 15 of the 20 
observed cycles should not be used. The purpose 
of this paper is not to attempt to evaluate the 
above but to present the information for the 
reader’s consideration. At any rate the reader 
may not find it difficult to believe that the use 
of the 19th cycle solar proton flux data may yield 
a pessimistic estimate of the solar proton hazard. 
With the above background, we can at least study 
the proton events characterized by the 19th cycle 
and infer proton events for future cycles similar 
to the 19th cycle. 

Next we consider the observation that during 
the 19th cycle, the large solar proton events did 
occur on the whole around the most active years of 
the 11-year cycle, however not necessarily in 
proportion to sun spot count. It is generally 
assumed that there were only about 6 years of 
observed large solar proton events. Thus the 
following analysis is based on the so-called 
active 6 years (300 weeks) of the cycle. The 
second assumption is that it is fair and logical 
to lump the actual dose rates over an active week 
into units of total dose per week (see Introduc- 
tion), This may include up to 3 proton events in 
a given week. Thus, instead of dealing with fun- 
damental units of proton events, we propose to 
deal with observed weekly dose rates behind 
various aluminum thicknesses. The data available 
from various sources gave a total of 24 proton 
events worth considering. (The 30 or more events 
discarded gave less than 3 rads total behind 5 g/ 


cm^ 0 f aluminum) . The results of grouping this 
data into weekly time periods gave a total of 18 
weeks with the frequencies of 13 weeks having 1 
event, 4 weeks having 2 events, and 1 week having 
3 events. Most of this data is recorded in NASA 
XND 4404 [3] on pages 16-23. The only major 
revision in reference 3 is in the use of the 
spectrum of A. J. Masleyt^] for the November 12, 
1960, flare. 

One note of explanation should be made re- 
garding the dose rates which are used in this 
study. They include a correction for secondary 
particles and thus are the sum of the primary 
proton dose and the secondary particles (neutrons 
and protons). This total dose at 20 g/cm 2 is 
about 20% and at 10 g/cm 2 about 10% above the 
primary proton dose rates. 

Now we come to a very treacherous part of the 
analysis. How do we choose a rational and suffi- 
cient model for an event dose-week for solar 
protons? Table 1 is presented as a summary of 
the weekly solar proton doses of the 19th solar 
cycle which the writer will consider. The absorbed 
doses in rads-tissue were calculated at the center 
point of a spherical shell of aluminum with the 
shell thickness being designated as the shield. 


TABLE 1: WEEKLY DOSES 07 THE 19TH SOLAR CYCLE 


Sample 

Number 



Number 

Shield 

6-Year Dose) j 


Week 

Year 

of 

Flares 

5(906)® 

101251) 

20(73) 

Remarks 

1 

Nov t2 - 18 

1960 

2 

456 

83. 4 

17. 0 

Category 1 - Large Weekly Doses 

2 

July 10 - 16 

1959 

3 

216 

75. 0 

21. 2 

0.01 

3 

Feb. 23 

1956 

1 

92 

50.0 

24.8 



Summed Dose 


(6) 

764 

208. 4 

63.0 

Approximately 85 Percent of 6-Year Total 


Percent of To- 
tal 6-Year Dose 



94.3 

S3 

86. 3 


1 

May to 

1959 

1 

59. 3 

18.3 

4. 4 

Category 2 - Medium Weekly Doses 

2 

July 18-20 

1961 

2 

22. 0 

8. 2 

2.4 

X 2 0.02 

3 

March 23 

1959 

1 

10.9 

2. 5 

0. 4 


4 

July 7 

1959 

1 

10.5 

2. 3 

n. 4 


5 

Jan. 20 

1957 

‘ 

8.3 

1. 8 

0. 3 


6 

Oct. 20 

1957 

1 

4. 1 

1. 8 

0. 7 



Summed Dose 


17) 

115. 1 

34.9 

8. 6 

Approximately 13 Percent of (i-Year Total 


Percent of To- 
tal (i-YearDose 



12.7 

13. 9 

11. 8 


t 

Aug. 22-26 

1959 

2 

5.9 

1.0 

0.2 

Category 3 - Small Meekly Doses 

2 

Aug. 29 

1957 

1 

4. 2 

0. 8 

0.1 

V 3 

3 

Nov. 20 

I960 

1 

3.6 

1. 5 

0. 05 


4 

Aug. 3 

1956 

1 

2. 2 

1. 0 

0, 4 


5 

Sept. 13 

1960 

‘ 

2.9 

1. 2 

0. 5 


a 

Aug. Ill 

1959 


1. 8 

0.4 

0. 1 


7 

July 3 

1957 

1 

1. 2 

0. 43 

0. 03 


« 

July 12 

1961 

1 

1.4 

0. 30 

0. 04 



April 2S - 29 

1*161) 

2 

0.77 

0. 27 

0 



Summed Dose 


Ml) 

23. 97 

6. 90 

1. 42 

Approximately 2 Percent of 0-Ycar Total 

L_ 

Percent of To- 
tal (i- Year Dost 




2.75 

1. 93 



The first number is the shield thickness in g/cm*-nluminum. 
the number in parenthesis is the O-year dose behind the shield. 
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The data of Table 1 have been grouped into 
three categories: ( 1 ) large weekly doses, ( 2 ) 
medium weekly doses, and (3) small weekly doses. 

Category (1), consisting of only 3 weeks, 
gave mory than 84% of the total 6 -year dose. 
Category (2), with 6 samples, gave about 13%, and 
Category (3), with 9 samples, gave less than 3% 
of the total dose. 

The grouping which is shown is certainly not 
unique, and the last 2 weeks of Category (2) per- 
haps should be in Category (3) , or perhaps the 
largest dose week of Category (2) should be in 
Category (1), The following work has attempted 
to provide for various combinations which the 
reader may wish to investigate by making the 
me thods of approach more important than the choice 
of a precise set of data. 

Since over 84% of the total solar proton 
dose is grouped into the 3 weeks of the first 
category, it would seem that one feasible solar 
proton dose model for an active week could be 
depicted by dividing the total 6 -year dose behind 
the various shields by 3. Since the most active 
phase of the 19th cycle is about 6 years, it seems 
sufficient to use 300 weeks as the basic period 
of major solar activity. With this assumption, 
which will be used throughout this study, it 
follows that the estimated weekly expectation of 
a dangerous large solar proton dose would be A 0 - 
.01 (event/week). 

In any case, the chance for an event of 
Category (2) or Category (3) is more likely to 
occur than one of Category (1) for a short ex- 
posure time. For this reason, a model will also 
be developed to reflect medium and small flare 
doses as well as large events. Thus, two solar 
proton dose models will be constructed and results 
compared in the following work. 

The next part of our approach is to choose 
the best composite 19th cycle solar proton flare 
data which is available and arrive at total tissue 
dose in rads from primary and secondary particles 
at the center of spherical shells with varying 
thicknesses. T. T. White, et. al^ have develop- 
ed a composite model (MSC model) for the total 
proton flux during the 19th cycle. This model 

gives a larger dose for energies less than about 
115 MeV than a model based on Webber's work[5]. 
Above 115 MeV, a composite flare model based on 
Webber's work (to be designated as the MSFC model) 
gives a larger dose. The 19th cycle composite 
6 -year proton spectrum which will be used in this 
work is the following: 

(MSC) I. J(>p) = 5.28 x 10 11 e “ p/73 

O0 < E < 115 MeV) (1) 

(MSFC) II. J(>p) = 1.14 x 10 11 e“ p/ *°° 

(E > 115 MeV) , 

where p is the rigidity units of MV and J(>p) is 
the integral spectra (proton/cm^) with energies 
above p. For protons, the relati onship betwe en 
p(MV) and E(MeV) is given by p = \/e^ + 1876E 
(ref. 3, page 15). 

Using the spectrum above, the best estimate 
of the 6 -year (300 weeks) total dose (primaries 
and secondaries) is given in Fig. 1. From Fig. 1 
and Table 1, the solar proton one-week-dose event 
models are constructed (Fig. 1 gives the magnitude, 
and Table 1 gives the fractions). Model I will be 
the dangerous solar event model which is represen- 
ted by 3 large events over a period of 300 weeks. 



The expected number of events in 1 week is given 
by *o = *01 (event /week) . Model II will depict 
the possibility of a large, medium, or small dose 
in a week where the percentages of the total dose 
in Table 2 are used to determine the relative 
size of a dose week in the three groups. For the 
three categories of Table 1, the values Ai = .01, 

X.2 = .02, and A 3 ** .03 are the mutually independent 
weekly expectations of large, medium, and small 
doses, respectively. Table 2 summarizes the dose 
models. It should be clear that the writer has 
presented in Table 2 only the primitive elements 
of a probability model. Thus the table gives the 
small sample estimates of the weekly expectations 
(A) and the consequences (rads/week) if an event 
occurs according to Model I or II. It should also 
be made clear that the values of A are valid only 
for the 300 most active weeks of a solar cycle 
similar to the 19th. The major unknown factor is 
an estimate of the chance of obtaining a solar 
cycle that would give total proton doses as large 
or larger than the 19th cycle. 

The values of A should be smaller for the so- 
called quiet sun (260 weeks) or else the magnitude 
of an event reduced to one of Category (3). Per- 
haps one might assume that A = .02 (one event per 
year) during the 260 weeks of the quiet sun, but an 
event should be depicted by the D 3 column (small 
event) of Model II. 

As a final conclusion to this section, it is 
of interest to recall that Fig. 1 represents the 
total dose versus spherical shell thickness for the 
19th cycle. If for the data of Fig. 1 as described 
above, an effort is made to correct for the self 
shielding of an astronaut and the complex geometry 
of a spacecraft which has aluminum walls of 13.5 
g/cm^ (assume that a 5 cm depth dose corresponds 
to an additional aluminum thickness of 6.5 g/cm^) 
then the 5 cm depth dose is estimated by dividing 
the dose at 20 g/cm^ by a factor of 3. Thus one 
finds that the total dose from the 19th cycle was 
25 rads at 5 cm tissue depth. Of course, this 
value is for the solar proton dose and does not 
account for the galactic cosmic ray dose which may 
range from 5 to 12 rads per year behind 20 g/cm^ 
depending on how much dose is contributed by the 
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heavy cosmic ray particles (Z > 2) . The impli- 
cations of all the above are simply that since the 
19th solar cycle was an example of a very active 
sun, and if one believes that the high energy pro- 
ton fluxes received during this cycle are associa- 
ted with this activity, he must be cautious in 
drawing certain conclusions. For example, the 
same astronaut above may have received from 30 to 
70 rads from galactic cosmic rays during a 6-year 
trip whereas the above dose from solar protons at 
bone marrow depth was 25 rads (13.5 g/cm 2 - al. 
walls) . Here one assumes that the galactic cosmic 
ray spectrum is so energetic that geometry factors 
are negligible for dose reduction. Thus the 
cosmic ray dose component may determine the limi- 
ting dose factor for long duration space travel. 

The purpose of the foregoing is not to minimize the 
importance of solar proton events but to point out 
that it may be quite feasible to shield against 
solar protons but probably impractical to consi- 
der shielding against galactic cosmic rays, and 
this very fact may be more important in determining 
man’s exposure time to space outside the earth’s 
magnetic field than the solar protons. 


probability that more than one event occurs is 
given by P n>1 (At) = ^(At). 

These conditions are the basic assumptions of 
the Poisson process (Fellert?], pages 400-402). It 
should be clear that we are stating that for a 
small time interval At, the chance for one event 
is approximately XAt, and the chance for more than 
one event is very small compared to XAt. Since t 
is for a relative time scale, it is not contradic- 
tory that At of 1 week can be small on our time 
scale. The above conditions lead to k system of 
differential equations for P n (t). They are: 

dP (t)/dt = -XP (t) + XP n > 1 , (4) 

n n n-± — 

and 

dP (t )/dt = -XP (t), n = 0 . (5) 

o o 

From Equation (5) and P q ( 0) = 1> we g et: 

P (t) * e ^ . (6) 

o 


A PROBABILITY MODEL IS DERIVED 

In this section, we shall address ourselves 
to the problem of choosing and using a probability 
density function in order to arrive at the proba- 
bility of getting ”x" events in ”t” weeks and the 
associated problem of probability or percentile 
levels. The important question of establishing 
confidence intervals for the basic statistical 
parameters which are obtained from a small sample 
will be undertaken in the next section. In this 
section, we will be making the naiye assumption 
that the basic statistics or population parameters 
(mean and variance) are well known, either by 
experience or apriori knowledge. 

In order to derive the probability model, 
which seems to be the most natural outgrowth of 
the solar proton dose week, the process will be 
described in terms of probabilities P n (t) that 
exactly n events occur during a time interval t 
(weeks in our case). Thus, P 0 (t) is the probabi- 
lity of no event in the interval t and 1 - P Q (t) 
is the probability of one or more events. Next we 
define X to be the mean or expectation of an 
event for a unit time interval. That is, 

£ = no. of events 

total weeks of observation 


is a statistical estimate of X. More precisely, 

X is a constant which determines the density of 
points on the t axis. Thus for a small interval 
of time At the probability of one or more events 
is given by: 

1 - P o (At) = XAt + £ (At) , (3) 

where e(At) is an infinitesimal and -small compared 
to XAt such that 


lim £ (At) 

At+o At 


(3a) 


Now we make the following postulate: 

Whatever the number of events during (o,t) the 
probability during (t, t + At) that one event 
occurs is given by P]_(At) = XAt + £ Q (At), and the 


Using Equation (6) and Pi(0) - 0, Equation (4) can 
be solved for Pi(t) * Xte“^ t . Using the fact that 
P n (0) = 0, (n>0) , Equation (4) becomes a recursion 
equation and successive values of P n (t) can be 
found. The resulting solutions give the terms of 
the Poisson distribution: 


P_(t) 


-Xt 


(Xt) n 


(7) 


If the reader accepts the postulate following 
Equation (3), the Poisson distribution is the 
natural outcome. The foregoing arguments have 
been presented in order to minimize the illusion 
that the author has pulled a distribution function 
out of the sky. For a more rigorous treatment, 
reference 7 is recommended. 

There are many interesting uses of the 
Poisson distribution in addition to the occurrence 
of rare events in a continuum of time. The 
distribution is used to approximate the binomial 
distribution for the case of rare events (p < .05). 
The word rare means individually rare. In a large 
population, several such events may occur, but the 
probability of occurrence of each individual event 
is small. For example, the number pf people 
killed by horses in 1969. An important feature 
of the Poisson distribution is that for large 
values of the mean (Xt » 20), the distribution 
approaches the normal (or Gaussian) distribution. 
There are many applications of the Poisson distri- 
bution given in any standard text on probability 
and statistics. The most common include such 
studies as the number born blind in a large city, 
radioactive disintegration, bacteria on plates, 
telephone traffic, etc. The remainder of this 
section will be devoted to a discussion of the 
Poisson distribution and how to apply it to our 
class of problems. 

If a discrete variable x has a Poisson dis- 
tribution, then the probability that [X - x] is 
given by the expression: 

-m x 

P[X = x] = — - , x = 0, 1, 2, (8) 


and P [ X = 0] = e m , since 0! - 1. 
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The mean (or expected value of the variable x) is 
given by y = m, the variance a 2 is given by a 2 =» m, 
and the standard deviation is simply given by a = 
VnT. The summation equation, 

CO 

E e~ m m X /x! * 1 , (8a) 

x~0 

is satisfied. Since we are primarily interested 
in the time-dependent form of the Poisson, m = At 
becomes the mean. Even though the Poisson is dis- 
crete in the variable x, it is continuous in t. 

For our purpose, we shall write the Poisson dis- 
tribution as: 

P[X - X ] = -■ ^, (Xt)X , x = 0, 1, 2, .... (9) 

where t = number of weeks, x - number of events 
in t weeks, A = mean number of events per week 
with (At) becoming the expected or mean number of 
events in t weeks, and P[X * x] is the probability 
of exactly x events in t weeks. 

If one wished to find a value of x which 
would not be exceeded say at least 99% of the 
time (1% chance at most of exceeding x) for a 
given At, the accumulative distribution function 
P[X < x*j is utilized. The value of x which we 
seek is the smallest x* which satisfies 

x* -At X x 

Pr [X < x*] = I - ^ P > (10) 

X=0 x- 

where P = .99 for the case in question. 

Since we are dealing with a discrete distri- 
bution function, the inequality is necessary, and 
one usually obtains a value of x* which corresponds 
to a value slightly above P (.99). This discre- 
pancy is usually circumvented by the proper use 
of the words "at least" and "at most." 

Next, we will examine how to apply the above 
information to our problem of predicting the solar 
proton dose that would be seen on a space mission 
at various percentile levels. For simplicity, let 
us assume that if an event occurs, it is depicted 
by the so-called "large" event model or Model I 
of Table 2. Using this model, the sample statis- 
tical estimate of A 0 « .01 event/week, and for 
convenience, assume the duration of exposure to 
large solar events is 100 weeks, then the value of 
At = 1.0 for the Poisson distribution. Now assume 
that we wish to find the dose levels that would be 
exceeded 27%, 2%, and 0.1% of the time. Applying 
Equation (10) with P taking on the values 0.73, 

0. 98, and 0.999, we see that the values of x are 

1, 3, and 5 events. So to find the dose at the 
respective levels, we multiply the doses behind 
the various shield thicknesses of Model I by the 
values of x above. The results are plotted in 
Fig. 2. Note that the 73.6 % level corresponds in 
our model to the mean or expected value of the 
Poisson distribution function. The percentile 
level at the expected value is found only when At 
is an integer. It becomes smaller approaching 50% 
as At increases to large values. For example, 
when At - 20, the percentile at x* - 20 is 55.9%. 

A similar computation can be carried out 
using event Model II. The values of the Poisson 
mean become respectively A^t - 1.0, \ 2 t = 2.0, and 
A 3 t = 3.0 where t 55 100 weeks t Tables are con- 
structed, and the following results are found at 
the expected value, 98% and 99.9% level. For the 


TABLE 2: MODELS OF A ONE-WEEK DOSE EVENT FOR THE ACTIVE 300 WEEKS OF THE 19TH 

SOLAR CYCLE 


MODEL I MODEL II 



N - 3* 
100X** 

X 0 - .01+ 

X(g/cm 2 ) 

D ) 

o Weeks/ 

2 

1420 

5 

340 

10 

93 

15 

43 

20 

25 

30 

11 

40 

6 

60 

2 



N(Large)-3 

84% 

X 1 - .01 

N(Med.)-6 

13% 

X 1 - .02 

N(Small)-9 

3% 

X 3 - .03 

3 

E NiDi 

1-1 

X(g/ctn 2 ) 

1 ^ Week ) 

n ( Raids \ 
°2 (v^ik) 

D, 

3 ( Week ) 

Total Dose 
300 Weeks 

2 

1195 

95 

15 

4290 

5 

286 

24 

3.4 

1032.6 

10 

78 

7.5 

0.9 

281.1 

15 

36 

3.5 

0.43 

132.87 

20 

21 

2.0 

0.25 

77.25 

30 

9.3 

0.80 

0.11 

33.69 

40 

5.03 

0.45 

0.06 

18.33 

60 

1.70 

0.15 

0.02 

6.18 


* The N indicates the number of weeks used in dose category. 

** The X gives the percent of total dose in the category, 
t The "X" gives the expectations (probability) of category for one week. 
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large events, x = 1, 3, and 5; for the medium 
events, x ■ 2, 5, and 8; and for the small events, 
x = 3, 7, and 10. The percentile values at the 
expected or mean value levels (x* = At) are 73.6%, 
67.7%, and 64.7%, respectively. In order to reduce 
the necessity of repetition, this set of levels 
could be denoted as approximately 70% level. The 
comparison of results from these totals to Fig. 

2 [8] indicate very good agreement. A useful rule 
would be: Empirical Rule - If at a given probabi- 

lity level there is at least one large event, then 
Model I is sufficient to describe the radiation 
hazard. If no large event is found at the given 
probability level, then Model II should be used. 

At this point, it might be useful to depict 
several curves similar to Fig. 2 for various 
mission times t, with A = .01 event/week, at 
given probability or percentile levels. In order 
to do this, arrays could be constructed for 
various values of At, and the number of events, 
x, necessary to satisfy the inequality of Equa- 
tion (10) could be found for various values of 
P < 1. 

Rather than provide a multitude of similar 
graphs, Table 3 will allow the reader to find the 
number of events for a range of anticipated 
values of At at 9 different percentile levels 
from 30 to 99.9%. Table 3 was constructed by 
choosing a probability level, P, the number of 
events, N, and then finding the value of m(-At) 
that satisfied the following equation: 

N -m x 

Pr(x £ N) = Z — — = P, (m = At) . (11) 
x=0 


TABLE 3: VALUES OF Xt TO SATISFY Pr(X j H) ■ ? 



0. 001 

0. 005 

0.010 

0.025 

0.050 

0. 100 

0. 15 

0. 25 

0.50 


99.9 

99. 5 

99.0 

97. 5 

95. 0 

90. 0 

85.0 

75.0 

50, 0 


Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

0 

0. 0010 

0. 0050 

0.010 

0.025 

0. 051 

0. 105 

0. 16 

0. 29 

0. 69 

1 

0. 045 

0. 10 

0. 15 

0. 24 

0. 36 

0. 53 

0. 68 

0. 96 

1.68 

2 

0, 19 

0. 34 

0. 44 

0. 62 

0. 82 

1. 10 

1.33 

1.73 

2. 67 

3 

0. 43 

0. G7 

0.82 

1.09 

1. 37 

1.74 

2.04 

2.54 

3. 67 

4 

0.74 

1.08 

1. 28 

1. 62 

1.97 

2. 43 

2. 79 

3. 37 

4. 67 

5 

1. 11 

1. 54 

1.79 

2. 20 

2.61 

3. 15 

3.56 

4. 22 

5. 07 

6 

1.52 

2. 04 

2. 33 

2. 81 

3. 29 

3. 89 

4. 35 

5. 08 

6. 67 

7 

1.97 

2.57 

2.91 

3. 45 

3. 96 

4. 66 

5. 15 

5. 96 

7.67 

8 

2. 45 

3. 13 

3. 51 

4. 12 

4.70 

5. 43 

5.97 

6. 84 

8. 67 

9 

2. 96 

3.72 

4. 13 

4. 80 

5. 43 

6. 22 

6. 80 

7.73 

9. 67 

to 

3. 49 

4. 32 

4.77 

S. 49 

6. 17 

7. 02 

7. 64 

8.62 

10. 67 

11 

4.04 

4. 94 

5. 43 

6. 20 

6, 92 

7. 83 

8. 48 

9. 52 

11.67 

12 

4.61 

5. 58 

6. 10 

6. 92 

7, 69 

8. 65 

9. 34 

10. 42 

12. 67 

13 

5. 20 

6. 23 

6. 78 

7. 65 

8. 46 

9. 47 

10. 19 

11.33 

13. 67 

14 

5. 79 

6.89 

7. 48 

8. 40 

9. 25 

10. 30 

11.06 

12.24 

14. 67 

15 

6. 41 

7.57 

8. 18 

9. 15 

10.04 

11.14 

11.92 

13. 15 

15. 67 

16 

7. 03 

8. 25 

8. 89 

9. 90 

10. 63 

11. 98 

12.79 

14.07 

16. 67 

17 

7. 66 

8. 94 

9. 62 

10.67 

11.63 

12.82 

13.67 

14. 99 

17. 67 

IB 

8. 31 

9. 64 

10. 35 

11. 44 

12. 44 

13. 67 

14.55 

15.91 

18. 67 

19 

8. 96 

10. 35 

11.08 

12. 22 

13. 25 

14. 53 

15.43 

16. 83 

19.67 

20 

9. 62 

11.07 

11.83 

13. 00 

14.07 

15. 38 

16. 31 

17. 75 

20.67 

21 

10. 29 

11.79 

12. 57 

13.79 

14. 89 

16. 24 

17. 20 

18.68 

21.87 

22 

10. 96 

12.52 

13. 33 

14. 58 

15.72 

17. 11 

18. 09 

19,61 

22. G7 

23 

11. 65 

13. 26 

14. 09 

15.38 

16. 55 

17. 97 

18. 98 

20. 54 

23. 07 

24 

12. 34 

14.00 

14.85 

16.18 , 

17. 38 

18. 84 

19. 88 

21.47 

24. 67 

25 

13. 03 

14.74 

15. 62 

16.98 

16. 22 

19.72 

20. 77 

22. 40 

25. 67 

26 

13.73 

15. 49 

16.40 

17.79 

19. 06 

20. 59 

21.67 

23. 34 

26.67 

27 

14. 44 

16. 25 

17. 17 

16.61 

19.90 

21.47 

22.57 

24. 27 

27. 67 

28 

15. 15 

17. 00 

17.96 

19. 42 

20.75 

22. 35 

23. 46 

25. 21 

28. 67 

29 

15. 87 

17.77 

18.74 

20. 24 

21. 59 

23.23 

24. 38 

26. 15 

29. 67 

30 

16. 59 

18. 53 

19. 53 

21. 06 

22. 44 

24. 11 

25. 23 

27. 09 

30. 67 

31 

17.32 

19. 30 

20. 32 

21.89 

23. 30 

25. 00 

26. 19 

26,03 

31.67 

32 

18. 05 

20. 08 

21. 12 

22. 72 

24. 15 

25. 89 

27. 10 

28.97 

32. 67 

33 

16. 78 

20. 86 

21, 92 

23.55 

25.01 

26.77 

28. 01 

29.91 

33. 67 

34 

19.52 

21. 64 

22.72 

24. 38 

25. 87 

27. 66 

28.92 

30. 85 

34. 67 

35 

20. 26 

22. 42 

23.53 

25. 21 

26. 73 

26. 56 

29. 83 

31.79 

35.67 

3G 

21.00 

23. 21 

24. 33 

26.05 

27. 69 

29. 45 

30. 75 

32.74 

36. 67 

37 

21.75 

24. 00 

25. 14 

26. 89 

28. 46 

30. 34 

31. 66 

33.68 

37. 67 

38 

22.51 

24. 79 

25. 96 

27. 73 

29. 33 

31. 24 

32. 59 

34.63 

38. 67 

39 

23. 26 

25. 59 

26.77 

28. 58 

30. 20 

32. 14 

33.50 

35.57 

39. 67 

40 

24. 02 

26. 38 

27. 59 

29. 42 

31.07 

33.04 

34.42 

36. 52 

40. 67 


a. P 1 - '* 

b. H At lies between two entries, use N corresponding to the largest entry. Thus, If 
At 0. 5, use N - 4 at the 3S. 3-percent level. If At is less than or equal to the lop 
entry, use N 0. Thus, if At 0.02, N - 0 at the 97.5-percent level, but N 1 at 
the 39. 0-percent level. 


The values of m - At, found in this manner, is 
the expected value or mean of the Poisson distri- 
bution which has a percentile level of exactly 
P x 100% for N or less events. 

The top of Table 3 is headed with a row of 
values labeled a which denotes the probability 
of more than N events or simple a - 1 - P. The 
probability statements using a will be made as 
follows: the probability is no greater than a 
that more than N events will be observed in t 
weeks. The use of a probabilities will be de- 
rived in the next section. The range of values 
in Table 3 should provide for the refinements 
needed in the following work. 

CONFIDENCE INTERVALS ARE ESTABLISHED 

The results of the previous section depended 
strongly on the choice of the parameter A or the 
probability of having a large event week. In 
this section, we will establish the confidence 
we can place on the value of A as calculated 
from observed sample data. One of the virtues of 
the Poisson distribution is that the value of the 
mean (m « At) completely determines the distri- 
bution function; i.e., » m, whereas the mean 

and variance (a^) is needed for most distribution 
functions, and they are not simply related. Thus 
in dealing with small sampling statistics from a 
normal distribution, one needs to establish con- 
fidence intervals for both the mean and the 
variance. 

In order to more clearly explain the intent 
of this section, an example will be given. The 
Poisson law arises very often in certain biologi- 
cal problems such as organisms distributed at 
random over the bottom of a lake. The number 
of such organisms found in a series of trial 
dredgings from separate small areas of the same 
size will follow this law. Statisticians cal- 
culate boundaries of possible outcomes from a 
given small sample, and these values are called 
confidence limits at a certain probability level 
for the assumed distribution function. Now, if a 
biologist counted 21 organisms from one of his 
dredgings, he could assert that he is 95% confi- 
dent that the mean or expected value lies between 
13 and 32 organisms per unit area assuming a 
Poisson distribution. Thus with only one sample 
and assumption of a Poisson process, it is possible 
to set upper and lower bounds on possible outcomes 
at a given probability level. From the above 
example, we can assert that if many dredgings are 
made, we expect only 5% will contain a number of 
organisms outside the predicted range. 

Returning to the basic problem, we wish to 
establish the confidence interval on m, hence to 
find the probability 

Pr(m" £ m £ m’) £ (1 - 23) 100% , (12) 

where 1 - P = 3. 

Now if P 5 .99 and 3 = .01, then 

Pr(m" £ m £ m f ) > 98%, (12a) 
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and we are at least 98% sure that m lies between 
m 1 and m", where m has only a 1% chance of being 
greater than m 1 . We can now write the general 
set of equations which will determine the value 
of m* and m": 



where x Q is observed number of events. The form 
of the above equation is an outgrowth of using a 
discrete distribution function and follows the 
common practice in textbooks on statistics. 

Solving Equations (13) and (15) for m f and m" 
by use of Newton’s method, the upper and lower 
bounds on m are found for "x 0 " observed events. 
Table 4 summarizes these results for 8 different 
values of 3 or P. Fig. 3 depicts a typical result 
taken from Table 4. 


TABLE 4: UPPER AND LOWER BOUNDS FOR THE MEAN 
AT PROBABILITY OF P = ( 1 - » ) 100% 



Even though we have obtained confidence 
bounds on the estimate of m * At, the major 
interest is concerned with the upper bounds that 
an investigator should use when only x Q events 
have been observed over a "t" week time interval. 
Following this concept and the simple model of 
a large dose- event -week, we recall that only 
3 such event weeks occurred over the 300 active 
weeks of the 19th solar cycle. Now if one wishes 
to find the upper bounds on this observation, he 
may use Table 4 in conjunction with Table 3 and 
arrive at the 100PZ confidence level that the pro- 
bability is no greater than some small value a 
that N^’ events will occur. Thus Table 5 is 



FIGURE 3: 902 CONFIDENCE LIMITS FOR POISSON 


TABLE 5: THE NUMBER OF EVENTS N' THAT WILL BE EXCEEDED 

a 


AT A PROBABILITY NO GREATER THAN a FOR THE 100-PERCENT 
UPPER BOUND VALUE OF THE POISSON MEAN, m‘, 

WHEN X - 3 OBSERVED EVENTS 


m ' 

v 

100P\ 

Percent 

0.001 

0.005 

0.010 

0.025 

0.050 

). 100 

0. 150 

0. 250 

0.500 

13. 06 

99.9 

25 

23 

22 

20 

19 

18 

17 

15 

13 

10.98 

99.5 

22 

20 

19 

18 

17 

15 

14 

13 

11 

10.05 

99.0 

21 

19 

18 

17 

15 

14 

13 

12 

10 

8. 77 

97. 5 

19 

17 

16 

15 

14 

13 

12 

11 

9 

7. 75 

95.0 

18 

16 

15 

14 

13 a 

11 

11 

9 

8 

6.68 

90.0 

16 

14 

13 

12 

11 

10 

9 

8 

7 

6.01 

85.0 

15 

13 

12 

11 

10 

9 

8 

7 

6 

5.11 

75.0 

13 

12 

11 

10 

! 9 

8 

7 

6 

5 

3.67 

50.0 

11 

9 

9 

8 

7 

6 

6 

5 

3 


a. Example: One is 95 percent confident that the true mean does not exceed 
7. 75 (X,j=3) , and using this value of the mean, one is certain that the 
probability of seeing more than 13 events during a 300 week active period 
is no greater than 0. 050. 

TABLE 6: THE NUMBER OF EVENTS N' THAT WILL BE EXCEEDED AT A 

Of 


PROBABILITY NO GREATER THAN a FOR THE 100P-PERCENT UPPER 
BOUND VALUE OF THE POISSON MEAN, m\ WHEN X - 0 
OBSERVED EVENTS 0 



V 

100P\ 










m 

Percent 

0. 001 

0.005 

0.010 

0.025 

0.050 

0. 100 

0. 150 

0.250 

0. 500 

6. 91 

99. 9 

16 

15 

14 

12 

11 

10 

10 

9 

7 

5.30 

99.5 

14 

12 

11 

10 

9 

8 

8 

7 

5 

4.61 

99.0 

12 

11 

10 

9 

8 

7 

7 

6 

4 

3.69 

97.5 

11 

9 

9 

8 

7 

6 

6 

5 

4 

3.00 

95.0 

10 

8 

8 

7 

6 

5 

5 

4 

3 

2.30 

90.0 

8 

7 

6 

6 

5 

4 

4 

3 

2 

1. 90 

85.0 

7 

6 

6 

5 

4 

4 

3 

3 

2 

1.39 

75.0 

6 

5 

5 

4 

4 

3 

3 

2 


0.69 

50.0 

4 

4 

3 

i 3 

2 

2 

2 

1 

0 
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constructed. 

Since the 20th solar cycle has produced no 
large events, it is of interest to ask if zero 
large events occur for a period of 300 active 
weeks, what are the possible upper bounds of m* 
and the number of events N a T that could be expected 
at some small probability. (This may be a nonsense 
question and will be discussed in the next 
section.) Table 6 is constructed with this in 
mind. The methods used are the same as for Table 
5. It is of interest to note that even though 
x Q - 0 events are observed that the 93% upper 
bound for the Poisson mean is m* * 3 which was the 
actual observed number of large events during the 
19th cycle. Thus it seems that for a reasonable 
level of confidence (95%) , the analyst would be 
justified if he used m* « 3 (X = 0.01 event/week) 
for large events even though the solar activity 
for a given cycle was considerably different than 
the 19th cycle. Perhaps after all is said, the 
only conclusion that can be drawn is that the ob- 
served 19th cycle dose events could be used for 
any near average solar cycle, and if a cycle is 
predicted to be similar to the 19th cycle, then 
the results of Table 5 (x Q = 3) should be seriously 
considered as a possible model. 

Finally, Table 7 is presented as a summary 
of values to use for X (events/week) at various 
confidence levels. For large events, it seems 
that the x Q = 0 column is probably reasonable to 
use if the solar cycle is not very active. If a 
cycle similar to the 19th is forecast, then the 
column under x Q = 3 is perf erred. 


TABLE 7: VALUES FOR X' - m* /300 (EVENTS/WEEK) FOR X 

o 

OBSERVED EVENTS AT 100P% UPPERBOUND CONFIDENCE LEVEL 


X 

100P%^^*s, 

0 

3 

6 

9 

99.9 

.0230 

.0435 

.0602 

.0755 

99.5 

.0177 

.0366 

.0522 

.0667 

99.0 

.0154 

.0335 

.0486 

.0626 

97.5 

.0123 

.0292 

.0435 

.0569 

95.0 

.0100 

.0258 

.0393 

.0524 

90.0 

.0070 

.0223 

.0351 

.0474 

85.0 

.0063 

.0200 

.0323 

.0442 

75.0 

.0046 

.0170 

.0285 

.0392 

50.0 

.0023 

.0122 

.0222 

.0322 

OBS* 

.0000 

.0100 

.0200 

.0300 


*0BS denotes the value of X /300 
o 


The value of X and the dose event model to 
use during the remaining 270 less active weeks of 
a solar cycle are not obtainable from the present 
analysis. However, until more data is available, 
the best one can do is use some essentially arbi- 
trary criteria. For example, if we believe that 
the chance of a large dose event week is definite- 
ly dependent on some minimal level of solar 
activity and this level is not approached during 
the quiet periods of the cycle, then it would be 
unreasonable to use the large event dose week even 
at a very low probability level such as the values 
shown under the x Q = 0 column of Table 7. Even if 
X is only .002 event /week for the Poisson distri- 
bution, one sees that the chance of getting two or 
more large events is approximately 6% for t - 
200 weeks. However, we cannot preclude the possi- 
bility of some smaller dose event occurring for the 
quiet period of a solar cycle. The present writer 
suggests using X = .02 event/week for the quiet 
period but recommends that the dose model for the 
small dose event (D3) of Model II (Table 2) be 
used. This is equivalent to expecting about one 
such event per year which is reasonably close to 
the actual observed number of events from October 
1961 through July 1966. 

In order to illustrate applications of the 
foregoing work, eight trip lengths will be con- 
sidered ranging from 13 to 260 weeks (5 years), 
during the most active 6 years of a solar cycle 
which is like the 19th cycle (xq =3). In order 
to simplify the possible combinations, four 
different probability levels for the upper bound 
mean (Xt) are used at four probability levels for 
the Poisson distribution. A mean (.01 t) corres- 
ponding to the actual observed (OBS) large event 
weeks of the 19th cycle is also given. This 
corresponds to the 95% level of the mean for x Q - 
0 events, and this value is recommended for the 
average type cycle. The results of these com- 
binations are shown in Table 8. The entries in 
this table give the number of large events that 
will be expected at a probability equal to or less 
than a corresponding to the Poisson mean (Xt) 
which we are 100P% sure will not be exceeded 
(see footnote on Table 5) . These tables are con- 
structed in a manner similar to Table 5, however, 
we start with Table 7 for values of X at various 
levels of percent confidence (100P) . 

Figures 4, 5, and 6 have been presented for 
values of a at .001, .01, and 0.1 for different 
upper bound values of the mean showing a range 
from the observed (approximately 50%) to the 99.9% 
confidence level of the mean (Xt) for a 78-week 
mission during the active weeks of a very active 
cycle. Comparisons are also shown in these 
figures to the work of other authors [2,9,10] w h 0 
have made similar computations but have used diff- 
erent models for prediction of the solar proton 
dose. An interesting aspect of the above compari- 
son is that at the 0.1% (100 xa) level, the 
present work is considerably lower even for the 
99.9% upper bound value of Xt. However, at the 
10% level, the reverse situation seems to be the 
case. The major difference in the methods of most 
of the other authors and the present work is that 
the size of our large event is fixed, but the 
number of events may be quite large, whereas in the 
other methods, the size of a single event may be 
extremely large (several times larger than the 
large event used in this work). These differences 
will be discussed in the last section. 
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A NEW PROBABILITY MODEL IS FOUND 


The foregoing section seems to provide some 
rationale for choosing a probability function for 
a given set of conditions: (a) the active 300 
weeks of a very active cycle such as the 19th, 

(b) the active 300 weeks of an average cycle, and 

(c) the quiet 270 weeks of any solar cycle. How- 
ever, if we examine Tables 5 through 8, a question 
presents itself. If we choose a very high confi- 
dence level for the Poisson mean and then deter- 
mine the number of events that would give a low 
probability a of getting a worse situation, then 
it seems that we are discussing a very improbable 
situation. That is, what is the chance of the 
mean being as high as the 99.9% confidence level, 
and if we use this mean the joint probability of 
seeing more than N a f events at a probability no 
greater than 0.001? For example, in Table 5, wc 
see that N a * * 25 events when P = .999 and a = 

.001. At first blush, one might suspect that the 
chance of both conditions occurring is on the 
order of 10“6. But care must be taken since we 
are dealing with cumulative distribution functions. 
Thus the product (1 - P)a does not correspond to a 
unique value of N a *. This can be ascertained by 
examining in Table 5 the values of (1 - P)a corres- 
ponding to N a * = 13. Now, we would like to ask 
the question; 

"If x c events are observed in a given time, t^, 
what is the probability of seeing more than N 
events in any observation time t?" 

With this knowledge a unique value of the proba- 
bility of seeing N events can be made for any 
period t _< T 0 . This question is also important 
because we propose to use the values of X calcu- 
lated for the total period of T n = 300 weeks and 
might be suspect of the true probability when one 
applies the value of X to a time period say of 
only 50 weeks. However, the reader surely agrees 
that the best estimate of the X corresponds to the 
total sample space of the 300 most active weeks 
during the 19th solar cycle. 

In order to provide an answer to the above, 
we must first ask what is the distribution of 
possible Poisson means if in a time t Q there are 
only x Q events observed. If we examine Equation 
(13), which gives the cumulative distribution of 
8 as a function of the upper bound values of m for 
a given x 0 , one sees that the derivative (-dft/dm) 
yields the desired upper bound probability density 
function for m; 

-m x o 

-d8/dm = f (m) - e y ■ . (16) 

o * 


This function is continuous in m 
called gamma distribution in the 
expected value of m, E(m) = x Q + 


O 1 (m) 


+ 1, and 


dm 


x ! 
o 


and is the so- 
variate m. The 
1; the variance. 


(17) 


To summarize the meaning of Equation (16) , 
we can state that for a given observation of x c 
events, the probability of the m being in the 
interval m to m + dm is given by the probability 
equation; 

-m x 0 , 

f (m)dtn = - — — — — . (18) 

o’ 


See Fig. 7 for illustration of Equation (18) with 
Xq - 3 . 



However, we wish to investigate the more 
subtle relationship that reflects the distribution 
of X's for a given observation of x Q events over 
a period of t Q weeks. Therefore, we make the 
change of variables denoted by m = Xt ? ; dm = t Q dX, 
and our probability density function in the 
variate X becomes: 

t x ° +1 X -t X 

f(x)dx = -Vr x ° e ° dX • (19) 

o* 

For a given value of X, the probability of seeing 
exactly x events in time "t" is given by the 
discrete Poisson distribution function; 

P(x) = Pr (x;Xt) = (xt ^, e , x=0,l,2 , . . . (20) 


where the probability of having a value of X in 
the interval X + dX is given by the density 
fuction of Equation (19) and t <_ t 0 (actually 
the case t > t D is equally valid) . 

Using the above definitions, the relation- 
ship we seek is given by the following: 

00 

Pr(x, t | x q , t Q ) = J Pr(x;Xt) . f(X) dX . (21) 

o 

Thus we are stating that for Equation (20) of the 
Poisson distribution, the probability of each 
possible \ (a spectrum of possible means) be 
folded into the equation and the results inte- 
grated over all possible values of X from zero to 
infinity. The result is a probability density 
function which is the expected value of the 
Poisson distribution over all possible means X. 
Thus , 

x +1 

t ° t X 00 x +x -(t +t)X 

Pr(x,t|x ,t )= ■ - 0 — j J X ° e ° dX. (22) 

* 1 o’ o x !x! J 
o o 


After integration, we obtain; 


x +1 

O X 


Pr(x,t|x o ,t Q ) 


t t r 

o 

x !x! 
o i 


(x^+x) ! 

( t +t ) x ° +x+1 
o 


(23) 


Now if we make the substitution 0 = t/t and sim- 
plify, the results are the discrete distribution 
function in the variate x: 


Pr(x|x o ,0) 


(x+x ) ! 
o 

x ! x ! 
o 


0 


x 


( 1 + 0 ) 


x+x 0 +l 


(24) 
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It can be shown that for the discrete distribution 
function above; , 

00 (x+x ) ! x 

1 — rrm - “ 1 

x=0 


x! x ! 
o 


(l+6) X+X o +1 


(25) 


and the mean or expected value of x is given by: 

x = E(x) = (x +1)0 . (26) 

o 


The proof of Equation (25) can be shown by resor- 
ting to hypergeometric functions, and the results 
of Equation (26) can be found by multiplying the 
summand of Equation (25) by x and after simplify- 
ing, the value (x Q +l)0 can be factored out leaving 
a sum that is the equivalent to that shown in 
Equation (25) • An illustration of the density 
function [Equation (24)] is shown in Fig. 8 for 
0=1 and 0 = .52, 


0 26 
0.24 
0.22 
0.20 
0.16 
0.16 
014 
012 
0.10 
0.06 











/ 

\ 





/ 

\ 

\ 




/ 

/ 


e “ 

V 


□ — X 0 -3. 0-1. TT-4 
O — X 0 -3. 0»62. V-2.08 

T THE DOTTED LINES ARE NOT VALID FO 
A DISCRETE DENSITY FUNCTION BUT A 




\ 

\ 


Se 




\ 

\ 


BE1 

rWEEN THE 

rwo cases ’ 

>HOWN 

_ 



/ 


\ 




/ 

/ 

_□ 

\ *1 

\ 

> 

k 

V 







/ 


% 

\ 

\ 

N 

l s 




/ 

/ 

/ 



N 

\ 



7 

/ 

/ 




\ 

\ 






















"'“'1 


0.02 

0 




1 



9 

7 1 



X EVENTS 

FLC'JFJL i: TilS SISCSETE FSCSABILITT DEHSITi Fu»CI Luii rrfij* ,a) 


This report is not intended to be a study of 
probability theory, and the ramifications of the 
above distribution function and its possible para- 
llels in other statistical work will not be pur- 
sued. 

Because of the format of Tables 5 and 8 and 
the manner in which we have previously made proba- 
bility statements, we wish to find the probability 
that more than N events are seen, given x Q and 0; 
or 

Pr(x>N|x o ,0) - 1 - Pr(x<N|x o ,G) , (27) 

where 6 = t/t Q> x 0 = number of events observed in 
time t Q (300 weeks), and t is the observation time 
during which N events are seen. 


Thus, the ultimate relation which we wish to 
answer our probability questions is given by 
Equation (27) or the obvious variations associated 
with it. In fact, if the foregoing is valid, we 
may dispense with the difficulties of choosing an 
upper bound value of the Poisson mean at a given 
level and then determining the probability that 
N events will be exceeded at some probability a as 
shown in Tables 5 and 8, where the true probability 
is actually not known. 

Using Equation (27), the number of combina- 
tions of 0, x Q , and N can readily get out of hand. 
For this reason, only the values of x Q observed 
events from 0 through 9 were used with different 
values of 0 from 13 weeks to 300 weeks (0 = 1). 

The values of N were extended to the point where 
Pr(x>N|xo,0) < 10“6. These tabulations are given 
in reference 8. 

It is very interesting to note that when 
x 0 * 0, 9 - 1, that the chance of seeing more than 
3 events is as high as .0625. This infers that 
even though no events were observed during a given 
300 week active cycle, we cannot be more than 
93.75% confident that 3 or less events could 
occur in a similar cycle. Also, we see from the 
same above assumption that we are 50% sure that we 
will see more than zero events. This seems to in- 
fer a dilemma bordering on the naive statement 
that if you know nothing about the probability of 
an event, you can only be 50% sure that it can’t 
happen; e.g. probability of life on Mars. This 
last statement seems to cast doubt on the useful- 
ness of the case when x Q = 0. 

However, it is of academic interest to inves- 
tigate further the case when x 0 = 0. For example. 
Equation (24) becomes: 

Pr (x| o,0) = 0 X /(1+6) X+1 (28) 

and for 0 > 0, we see that when x = 0, 

Pr(x = o|o,0) = 1/(1 + 6) , (29) 

and 

Pr(x > o|o,0) = 0/(1 + 0) . (30) 

The probability of seeing an event as 0 approaches 
zero becomes very small as one would suspect for a 
very rare event. Now as 0 increases to large 
values (0 » 1), the value of the probability 
approaches 1. This infers that if an event can 
happen at a given small probability, then it is 
almost certain that the event will occur after a 
sufficiently long period. Thus Equations (29) and 
(30) do not defy intuition in an ordinary sense 
but leave us with a rather insecure feeling since 
the number of actual observed events in time t Q is 
zero. However, if we examine the density function 
of possible values of X [Equation (19] when x 0 = 

0, a plausible probability density function is 
found; 

f(X) dX = t e _t< A dx . (31) 

Table 9 is presented as a survey of Equation 
(27) for Pr (x > N | 3 , 0) <^ £ at several values of £ 
and 0. 
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A point of interest concerns the values of 
the results in Table 8 as compared to those in 
Table 9. It should be clear that the results of 
Table 9 do not provide the same type of probabili- 
ty statements as Table 8 . In general, this 
writer feels that the results of Table 9 are more 
useful and to the point but others may disagree. 

In any case, one may ascertain which probability 
level of the parameter X is more credible when 
using the methods of the previous section and con- 
structing tables similar to Table 8 . For example, 
one sees that the value of X at the 90% level 
seems to yield event numbers (N) at probabilities 
a in Table 8 which seem to be comparable to the 
probabilities e in Table 9, 

Finally Table 10 is presented for the case 
when x Q = 0 even though we have cast doubt on the 

validity of the meaning of this rather extreme case. 
However, it does signify a sort of boundary con- 
dition for those periods when the probability of a 
large event is very small as perhaps exists for 
the active years of an average type solar cycle. 

REVIEW AND COMMENTS 

In the foregoing sections, the writer has 
attempted to convey a method of th'inking about the 
radiation hazards associated with solar proton 
events. The method of approach is felt to be more 
important than the actual results presented. The 
methods and models as developed utilize very little 
solar physics as such and consequently will be very 
unsatisfying to many physicists who have examined 
various aspects of the problem. There is no 
attempt to model or predict a solar proton spectrum, 
the time dependence of particle arrival, or angular 
distributions of the flux. 

As a brief summary, the total dose in rads- 
tissue from solar protons was tabulated for weekly 
time intervals , and the number of weeks which gave 
a dose above 25 rads behind 10 g/cm^ of aluminum 
for the active 6 years of the 19th cycle were 
called dangerous or large event weeks. The number 
of such event weeks was found to be only 3 weeks 
for the past 20 years. 


Even though the chance for smaller events is 
examined, it was found that for any reasonable high 
confidence level (95%) , the smaller events could be 
ignored. Consequently, we took the total particle 
flux for the 19th cycle and divided this spectrum 
by a factor of 3 and arrived at a single large 
event week. 

Using this spectrum, one can calculate the 
tissue dose in rads at the center of an aluminum 
spherical shell (Table 2). To correct for geo- 
metric effects and self shielding, this dose should 
be reduced by a factor of about 3. If the space 
mission is planned during the quiet period of a 
cycle, then the small event dose curve (D 3 ) of 
Model II (Table 2) may be used with X * .02 (j* - 1 
event /year) . 

To predict the probability of an event 
occurring, the Poisson distribution seemed to be 
the most logical choice. The third section was 
devoted to examining this conclusion and the 
methods of using this probability model. The 
fourth section was written as an effort to define 
the confidence one could use in employing the 
Poisson process and specifically arriving at con- 
fidence levels for the experimental or observed 
value of the mean of the Poisson distribution 
function. Several examples were given for dif- 
ferent mission lengths, and comparisons were made 
to other authors' results. 

Finally, the previous section was an extension 
of the Poisson process to incorporate the concept 
of small sample theory and arrive at the expected 
distribution function which answers the following 
question. If x Q events are observed in time t 0 , 
what is the probability of seeing x events in any 
observation time t? The results were represented 
by the discrete probability density function in 
the variate x: Pr(x|xo,0). Using the above 

function, extensive tables were tabulated in 
reference 8 . 

In the beginning of this work, the author in- 
tended to avoid commenting on the methods used by 
other investigators, but in order to explain the 
radical differences shown in Figs. 4, 5, 6 , the 
following comments are pertinent. One of the more 
common procedures used is to obtain the logarithms 
to the base ten of the solar proton flux above 30 
MeV for each of the events of the 19th cycle. This 
data is then plotted on normal probability graph 
paper, obtaining a distribution called the log 
normal distribution. in the variable x ■ log^o 

This is a true normal distribution in the 
variate x, and there is no virtue in examining the 
transformed distribution in the variable <J> which 
has a more complex representation. Also the values 
of the mean and standard deviation between the two 
distributions are not simply antilog related. For 
more details, see reference 8 . 

The present writer has assumed that the users 
of this distribution keep their statistics in the 
variable x which is certainly the simplest process. 

From reference 10, 60 events were taken that 
had fluxes measured above 30 MeV. 

The logarithm of these 60 entries were ordered 
from the smallest to the largest and the normalized 
cumulative sums plotted on log normal paper as 
shown in Fig. 9. The ungrouped data had a mean of 
7.39, and the standard deviation was 0.97 as shown 
at bottom of the figure. The straight line in 
Fig. 9 depicts the cumulative normal distribution 
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TABLE 9: NUMBER OF EVENTS N THAT SATISFY Pr(X>N|3,0) £ e 


WEEKS 

Nv e 
0 

.10 

.01 

.001 

.0001 

.00001 

13 

.0433 

1 

2 

3 

3 

4 

26 

.0866 

1 

2 

4 

5 

6 

52 

.1733 

2 

4 

5 

7 

8 

78 

.2600 

3 

5 

7 

8 

10 

104 

.3467 

3 

6 

8 

10 

12 

156 

.5200 

4 

8 

10 

13 

16 

208 

.6933 

6 

9 

13 

16 

19 

260 

.8667 

7 

ii 

15 

i 

19 

22 

300 

1.000 

9 

13 

17 

21 

25 


with mean of 7,39 and variance of 0.97. from 
examination of this fit, one might conclude that 
the log normal gives a reasonable representation 
of the data. Note that the flux for the November 
12, 1960, event is plotted as a square which falls 
exactly on our line of best fit. This is done 
since actually when one obtains the cumulative 
distribution from a discrete set of data, the last 
point has the cumulative probability of 1.00, but 
since this cannot theoretically occur, the best 
choice is to place this point on the best fit line. 
This point which is logio <t> = 9.9562 corresponds 
to the cumulative probability of 99.6% or in terms 
of the standard deviation, the November 12, 1960, 
event is at 2.63 standard deviations above the 
mean. If we wished to go to the 99.9% level, 
then it would be necessary to take 3.09 standard 
deviations above the mean or a log flux of 10.4016. 
Now if we find the antilog of this value, we see 
that the flux at the 99.9% level in the log normal 
distribution is 2.52 x 1010 proton/cm 2 (E > 30 MeV) 
which is a factor of about 2.8 times larger than 
the November 12, 1960, event. Hence, it would 
seem that a flux above 30 MeV which is 3 times 
larger than the November 12, 1960, event would 
have a probability of occurring which is less than 

0.001. 

From the above analysis, one may be convinced 
that a reasonable upper bound value for a single 
event is at most a factor of 3 larger than the 
November 12, 1960, event. For the above reason 
this writer feels that the results obtained by 
some investigators for the extreme probability 
tails must depict a smaller probability than the 
estimates given. They have possibly used joint 
probabilities of flux and rigidity parameters which 
may be a factor of 10 or more smaller than those 
indicated by the 0.1% probability tail shown in 
the various reports at my disposal. 

If the reader wishes to use the 99.9% pro- 
bability event in the present work, he may ignore 
self shielding and geometry factors, a factor of 3. 


TABLE 10: NUMBER OF EVENTS N THAT SATISFY Pr(X> N|O,0) <Z 


t 

WEEKS 

N \ e 

.1 

.01 

' 

.001 

.0001 

.00001 

13 

.0433 

0 

1 

2 

2 

3 

26 

.0866 

0 

1 

2 

3 

4 

52 

.1733 

1 

2 

3 

4 

6 

78 

.2600 

1 

2 

4 

5 

7 

104 

.3467 

1 

3 

5 

6 

8 

156 

.5200 

2 

4 

6 

8 

10 

208 

.6933 

2 

5 

7 

10 

12 

260 

.8667 

2 

5 

8 

11 

14 

300 

1.0000 1 

3 

6 

9 

13 

16 
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